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, Abstract 

• The behaviour of expanding cosmological models with collisionless 

iy~j | matter and a positive cosmological constant is analysed. It is shown 

. that under the assumption of plane or hyperbolic symmetry the area 

' radius goes to infinity, the spacetimes are future geodesically com- 

. plete, and the expansion becomes isotropic and exponential at late 

| times. This proves a form of the cosmic no hair theorem in this class 

— «. . of spacetimes. 

cr 

tl 1 1 Introduction 

op. 

m ^ ' The presence of a positive cosmological constant A can lead to exponential 

^ . expansion in cosmological models. This is the simplest mathematical de- 

scription of an inflationary universe. Under certain circumstances the de 
Sitter solution acts as a late time attractor for more general solutions of the 
Einstein equations with A > 0. This is sometimes known as the cosmic no 
hair theorem. Up to now there are unfortunately not many cases where this 
kind of statement has been proved rigorously for inhomogeneous spacetimes. 

A positive cosmological constant can be introduced in Newtonian cos- 
mology and this provides a simplified model for the general relativistic case. 
In 4 a form of the cosmic no hair theorem was proved in Newtonian cos- 
mology. A perfect fluid was used as a matter model and solutions were 
considered which evolve from initial data which are small but finite pertur- 
bations of homogeneous data. It was shown that if the homogeneous solution 
exists globally in the future the same is true of the inhomogeneous solution. 
Of course a global in time existence theorem is a prerequisite for a proof 
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of the cosmic no hair theorem. It was then shown that the inhomogeneous 
solutions have a behaviour at late times which is qualitatively similar to that 
of the homogeneous model. If p is the mean density and bp = p — p then 
bp/p converges as t — > oo. In Newtonian cosmology there is also a theo- 
rem about the late time asymptotics of models with a kinetic description of 
matter by the Vlasov equation and with vanishing cosmological constant |Sj. 
The boundedness of bp/p is also obtained in that case. Adding a positive 
cosmological constant to the problem considered in 8 would presumably 
simplify the analysis but this has not been attempted. 

In general relativity the problem of proving the cosmic no hair theorem 
is more difficult. In the spatially homogeneous case there is a general result 
of Wald on spacetimes with positive cosmological constant which does 
not depend on the details of the matter content but only on energy condi- 
tions. There is one example where future geodesic completeness has been 
proved for a class of inhomogeneous spacetimes with matter |7]. This con- 
cerns solutions of the Einstein- Vlasov system with hyperbolic symmetry and 
A = satisfying an additional inequality on the initial data. Under those 
assumptions geodesic completeness was proved but only limited information 
was obtained on the asymptotic behaviour at late times. In the following 
we will show that in the presence of a positive cosmological constant this 
result can be strengthened a lot. Future geodesic completeness is proved 
for all solutions with hyperbolic and plane symmetry. Moreover the asymp- 
totic behaviour is shown to closely resemble that of the de Sitter solution. 
It should be mentioned that in the case of the vacuum Einstein equations 
there is a proof of a form of the cosmic no hair theorem which does not 
require symmetry assumptions but does require a small data restriction jSJ. 

In this paper we study solutions of the Einstein equations with positive 
cosmological constant coupled to the Vlasov equation describing collisionless 
matter. Under the assumption of plane or hyperbolic symmetry we show 
that the solutions are future geodesically complete and we obtain a detailed 
description of their late time behaviour, which is similar to that of the 
de Sitter solution. The proof is built on the local existence theorem and 
continuation criterion previously proved by Tchapnda and Noutchegueme 

Let us recall the formulation of the Einstein- Vlasov system which gov- 
erns the time evolution of a self-gravitating collisionless gas in the context 
of general relativity; for the moment we do not assume any symmetry of 
the spacetime. All the particles in the gas are assumed to have equal rest 
mass, normalized to unity. The four-momentum of each particle is a future- 
pointing unit timelike vector so that the number density / of particles is a 
non-negative function supported on the mass shell 

PM := {g a( 3p a ^ = -1, p° > 0}, 

a submanifold of tangent bundle TM of the space-time manifold M with 
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metric g (the signature is — h ++). We use coordinates (t,x a ) with zero 
shift and corresponding canonical momenta p a ; Greek indices always run 
from to 3, and Latin ones from 1 to 3. On the mass shell PM the variable 
p° becomes a function of the remaining variables (t,x a ,p b ) : 



where p a = g a /3P , are the Christoffel symbols, \g\ denotes the determi- 
nant of the metric g, G a p the Einstein tensor, A the cosmological constant, 
and T a a is the energy-momentum tensor. 

In this paper we want to investigate the Einstein- Vlasov system with 
positive cosmological constant and plane or hyperbolic symmetry in the 
expanding direction. For the notion of spherical, plane or hyperbolic sym- 
metry, we refer to jHj. We write the system in areal coordinates, i.e., the 
coordinates are chosen such that R = t, where R is the area radius function 
on a surface of symmetry. 

The circumstances under which coordinates of this type exist are dis- 
cussed in [2] for the Einstein- Vlasov system with vanishing cosmological 
constant. It will now be shown that the analysis there can be extended to 
the situation under consideration here. Consider first the Einstein equations 
with a general matter model satisfying the dominant energy condition and 
A = 0. For plane symmetric spacetimes it follows from Proposition 3.1 of 
jl()j that the gradient of R is always timelike. The corresponding statment 
in the case of hyperbolic symmetry can be proved by the argument in Step 
1 in section 4 of |2j. When there is a cosmological constant we can consider 
it as a fictitious matter field with 'energy- momentum tensor' —Ag a g. This 
fictitious matter field satisfies the dominant energy condition. The same 
is true of the tensor which is the sum of the fictitious energy-momentum 
tensor with the energy-momentum tensor of real matter satisfying the dom- 
inant energy condition. It can be concluded from all this that the gradient 
of R is timelike for the Einstein- Vlasov system with positive cosmological 
constant and plane or hyperbolic symmetry. The remainder of Step 1 and 
Step 2 in section 4 of [T] and [2| can be extended to the case where a posi- 
tive cosmological constant is present using the same method of the fictitious 
energy-momentum tensor. The only property which is required in addition 
to the dominant energy condition is the inequality q < p—p and this is satis- 
fied by the fictitous energy-momentum tensor. From this point it is possible 



P° = V^F \A + 9abP a P b - 



The Einstein- Vlasov system now reads 



d t f + ^d x af - l r ^yy«v/ = o 
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to argue exactly as in and |2j to conclude that a solution of the Einstein- 
Vlasov system with positive cosmological constant and plane or hyperbolic 
symmetry contains a Cauchy surface where R is constant. Hence there is no 
loss of generality in restricting consideration to spacetimes evolving from a 
hypersurface of constant areal time. 

Since the gradient of R is everywhere timelike it must be either every- 
where future-pointing timelike or everywhere past-pointing timelike. We 
choose a time orientation such that the latter is the case. Then the expand- 
ing direction of the cosmological model corresponding to increasing area 
radius t. 

The metric takes the form 

ds 2 = _ e 2»{t,v) dt 2 + e 2X(t,r) dr 2 + ^2 + ^2 0d(p 2 ) (1 -Q 

Here t > 0, the functions A and \x are periodic in r with period 1 and 

{sin# if k = 1 
1 if k = 

sinh# if k = — 1 

For the case k = 1 the orbits of the symmetry action are two-dimensional 
spheres. In this spherically symmetric well as in the case A < 

0, the global result below is seen to be false, cf. so these cases will 
not be considered further. For the plane symmetric case k = the orbits 
of the symmetry action are flat tori, they are hyperbolic spaces for the 
hyperbolic symmetry k = — 1. The coordinates (8, ip) range in [0, 2tt] x [0, 2tt] 
or [0, oo[x [0, 2tt] for k = 0, —1 respectively. It has been shown in [£] and |2j 
that due to the symmetry / can be written as a function of 

t, r, w := e V and F := t 4 (p 2 f + t 4 sin^ 6(p 3 ) 2 , 

1. e. / = f(t,r,w,F). In these variables we have p° = e^y/l + w 2 + F/t 2 . 
After calculating the Vlasov equation in these variables, the non-trivial com- 
ponents of the Einstein tensor, and the energy-momentum tensor and de- 
noting by a dot or by prime the derivatives of the metric components with 
respect to t or r respectively, the complete Einstein- Vlasov system reads as 
follows : 

dtf + -=jL=^L==d r f - (Xw + e »- x n'^l + w 2 + F/t 2 )d w f = (1.2) 

y'l + W Z + F/t A 

e~ 2fl (2t\ + l) + k- At 2 = 8nt 2 p (1.3) 
e^(2tfi-l)-k + At 2 = 8nt 2 p (1.4) 



4 



fjl = -4TTte x+fl j 



(1.5) 



e- 2X (// + //(// " " e~ 2M (a + (A - A)(A + + A = 4^ (1.6) 



where 

•OO /•OO 



P(*. r )-=72 / VT+w* + F/1?f(t, r, w, F)dFdw = e^T^t, r) 

t i-oo JO 

(1.7) 

/-oo /-co w 2 
* 2 J-oo Jo y / l + W 2 + F/t 2 

(1.8) 

/-oo /"OO 

j(tir):=-t / wf(t,r,w,F)dFdw = -e x+ > 1 T 01 (t,r) (1.9) 
1 J-oo Jo 

roo poo p 2 

?(*> r ) := 74 / / 7l 5 ^=f /& r, F)dFdw = ^T 22 (t, r). 

t A J-oo Jo y / l + W 2 + F/t 2 t 2 

(1.10) 

We prescribe initial data at some time t = to > 0, 

° ° o 

f{t ,r,w,F) = f(r,w,F), \(t ,r) = A(r), fi(t ,r) = /j,(r). 

The organization of the paper is as follows. In the next section we show 
that the solution of the initial value problem corresponding to the Einstein- 
Vlasov system (|1.2jl - (|1.6JI exists for all t > to- I n section 3 we first prove that 
the spacetime obtained in section 2 is timelike and null geodesically complete 
towards the future; later on we determine the explicit leading behaviour of 
solution, compute the generalized Kasner exponents and prove that each of 
them tends to | as t tends to +oo. 

2 Global existence in the future 

In this section we follow the approach of [2j . We make use of the continuation 
criterion in the following local existence result : 

Proposition 2.1 Let f G C 1 ^ 2 x [0,oo[) with f(r + l,w,F) = f{r,w,F) 
for (r,w,F) el 2 x [0,oo[, / > 0, and 

o 

wq := sup{\w\\(r,w, F) G supp/} < oo 
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o 

Fq := sup{F| (r, w, F) E supp/} < oo 



Let A G C 1 ^), pL G C 2 (R) with A(r) = A(r + 1), £(r) = Ji(r + 1) for r£l, 
and 



ft) = -4^ e A+ ^(r) 



"17 



3 A+M 



OO /"OO 



oo JO 



wf(r, w, F)dFdw, reK. 



Then there exists a unique, right maximal, regular solution (/, A, fj) of hl.ty) - 

° ° o 

\l.b)) with (/, A, /x)(to) = (/>A, fi) on a time interval [to,T[ with T E]to,oo\. 
If 

sup{/i(t,r)|r 6K,t£ [t ,T[} < oo 

then T = oo. 

This is the content of Thms. 3.3 and 3.4 in jllj . For a regular solution all 
derivatives which appear in the system exist and are continuous by defini- 
tion, cf. [TT] . 

We now establish a series of estimates which will result in an upper 
bound on \x and will therefore prove that T = oo. Similar estimates were 
used in for the Einstein- Vlasov system with Gowdy symmetry and were 
generalized to the case of T 2 symmetry in Unless otherwise specified in 
what follows constants denoted by C will be positive, may depend on the 
initial data and on A and may change their value from line to line. 

Firstly, integration of (|1.4|) with respect to t and the fact that p is non- 
negative imply that 



to(e 



> 



-Mr) +k) , 87T /■* 2 . N , A 3 ^ 
— '- -k--jr jf s 2 p(s,r)ds + -(i 3 - t 3 ) 



, *€ [t 0) T[. 



C - kt + ft 3 



In this inequality, C does not depend on A. Next we claim that 



f e^ +X p{t,r)dr < Ct~ 3 , t G [t ,T[. 
Jo 



A lengthy computation shows that 



4- C e» +X p{t,r)dr = -- C 
dtj PV ; t J 

Now using (|2.1|) we have 



(2.1) 



(2.2) 



2 p + g_£i^(i + A ;e ^_AtV'') 



<ir. 
(2.3) 



1 + ke 2 » - At 2 e 2 ^ < 1 + 



kt - At 3 



C - |At 3 



C-kt + ^t 3 C-kt + ^t 3 ' 
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The right hand side of this inequality is neg ative if t > (^) 1/3 . In this case 
1 + ke 2 ^ - At 2 e 2 ^ < so that, using the fact that q > and p > 0, (|!Hf> 
implies that 



d_ 

dt 



e^ +x p(t,r)dr < 



2p- |(1 + ke 2 ^ - Ai 2 e 2 ^)J dr. (2.4) 



Setting C"(A) := ;f(3C — 2k), we have the estimate 

1 + ke ^ _ At 2 e^ < 1 + „ kt 7 M \ ,„ < C'(A)t- 2 - 2, 



C — kt + ft 3 



and combining this with Q2.4JI yields 
which multiplied by t 3 gives 



e fl+x p(t,r)dr 









t 3 [ e^ +x pdr 


< Ct" 3 




. •/ 







d_ 
dt 

By Gronwall's inequality, this implies that 

ir < Ct- 3 , 



[ e» +x p{t,r)dr 
Jo 



i.e. (|2~2j) for i > (§x) 1/3 . For t < (^) 1/3 , (|2~3j) implies the following, since 
<?>0: 



f e^ +x p(t,r)dr < -- 



dt 



o 



t 



dr 



o 



+ 



< 



1 f 1 (k — Ai 2 )e 2 ^+ A 



e^pdr, 



(p + p)(ir 



we have used the fact that p > p and k — At 2 < 0. By Gronwall's inequality, 
we obtain 



e^ +x p(t,r)dr < Cr 1 

< (Ct 2 + C)t~ 3 

< 



t 3 since 



t <(|f>"» 
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that is ((221) nolds for t < (It) 1/3 

as well. Using (|2.2|) and the equation 

fjf = —4irte fl+x j we find 

pi pi pr pi pi 

| ju(t, r) — / fi(t, a)da | =| / / /j!(t, r)drda \< j / \n'{t, r)\drda 

JO JO Jo Jo Jo 

<4irt [ e^ +x \j(t,T)\dT <4irt [ e^ +x p(t,r)dT 
Jo Jo 

<Ct-\ te [to,T[, re [0,1]. (2.5) 

Next we show that 

ef*(*.r)-A(t,r) < Ct -2^ f £ r £ [ 0> ^ ( 2 g) 

To see this observe that by (fO)l . lfO)l and l(2~Tj) 
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1 + fee 2 ' 1 



Ate 2 ' 4 



< 



1 k - At 2 
- + 



/«-A 



Using the fact that — /c + Ai 2 is the derivative of C — kt + — t 3 and integrating 
this inequality with respect to t yields 



e^ x < C 



C-kt + ft 3 



< cr 



i.e. SFty . 

We now estimate the average of \x over the interval [0, 1] which in com- 
bination with (|2.5j) will yield the desired upper bound on /j, : 



/ /J<(t,r)dr = / 
Jo Jo 



r)dr 
t 1 rl 



t /•! 



to JO 



fi(s, r)drds 



<C+ [ — I \e 2 ^(87rs 2 p + k- As 2 ) + ljdrds 
Jt 2s Jo 



<C+~ln(t/t ) + C j s- A ds-^ 



to 



-k + As 2 



2 J t<] C-ks + f s 3 



In 



C - ks + f s 3 



where we used (|2~T)t . (j2~2)) . (|2~B|) and the fact that p < p. With lf23|) this 
implies 

fi(t, r) < C(l + t" 2 + In t~ 2 ) < C, i € [t , T[, r £ [0, 1] (2.7) 
which by Proposition 12. II implies T = oo. Thus we have proven : 
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Theorem 2.2 For initial data as in Proposition \2.1\ the solution of the 
Einstein- Vlasov system with positive cosmological constant and plane or hy- 
perbolic symmetry, written in areal coordinates, exists for all t G [to,oo[ 
where t denotes the area radius of the surfaces of symmetry of the induced 
spacetime. The solution satisfies the estimates 12. 12. fij) and \2. 7j ). 



3 On future asymptotic behaviour 

In the first part of this section we prove that the spacetime obtained in 
Theorem 12.21 is timelike and null geodesically complete in the expanding di- 
rection. The analogue of this result was proved by Rein, cf. 0, in the case 
A = 0, k = — 1 but with initial data satisfying a certain size restriction, an 
additional assumption which we drop here due to the fact that A does not 
vanish. The proofs of the results obtained in the first two subsections are 
modelled on the approach of [7j. In this section we are interested in prov- 
ing statements about the asymptotic behaviour of solutions at late times. 
Therefore there is no loss of generality in prescribing data at some large time 
t = t > 0. 

Firstly we establish a bound on w along characteristics of the Vlasov 
equation. 

3.1 An estimate along characteristics 

Let 

o 

wq := sup{\w\\(r,w, F) G supp/} < oo, 



Fo : = sup{F|(r, w, F) G supp/} < oo. 
Except in the vacuum case we have wq > and Fo > 0. For t > to define 
P+(t) := max{0, max{w|(r, w, F) £ supp/(t)}}, 

P-(t) := mm{0,mm{w\(r,w, F) G supp/(t)}}. 
Fix e G]0, 1[. We claim that 

P+(t) < wo (£) , P_(t) > -w (J^j , i > to. (3.1) 
Assume that the estimate on P + were false for some t. Define 

ti := sup | t > t \P + (s) < wq f^j ,t < s < 1 1 
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so that t < ti < oo and P+(h) = w (j^J > 0. Choose a e]0, 1[. By 
continuity, there exists some ti > t\ such that the following holds : 

(l-a)P+(«) > 0, s e [ii,i 2 ]. 

If for some characteristic curve (r(s), w(s), F) in the support of /, that is 

o 

with (r(to),w(to), F) S supp/, and for some t E]ii,i2] the estimate 

(1 - a/2)P+(t) < w(t) < P+(t) (3.2) 

holds then 

(1 - a)P+(s) < w(s) < P+{s), s G [h,t]. (3.3) 

Note that the estimates on w from above hold by definition of P+ in any 
case. Let (r(s),w(s), F) be a characteristic in the support of / satisfying 
(|3.2|) for some t G]ii,i 2 ] and thus (|3.3|) on [ti,t\. Then on [ti,t], 

w = —e 2 ^ / ( w^Jl + w 2 + F/s 2 - w\] 1 + w 2 + F/s 2 J /dFefo) 

1 + fee 2 " A 2 

H w swe p 

2s 2 

< !l[! e 2M / P+(S) / F ° ^ 2 (l + ^ 2 + F/s 2 ) - ^ 2 (1 + u> 2 + F/s 2 ) 

w 



s Jo Jo 
1 + ke 2 ^ A 



y/l + w 2 + F/s 2 + Wl + w 2 + F/s 2 



H : swe 2>l 

2s 2 



< / / — -fdFdw H w swe Zfl 

s Jo Jo w 2s 2 

.2^/ m x „ °, „ e 2 ^ , 1 1 + A;e 2 ^ - As 2 e 2 ^ , „, 

< 4vr 2 F (l + F ) || / || 27^(0- + ^ «> ( 3 - 4 ) 

o 1 1 _|_ t-pS/x _ A 2 2/i 

< 4vr 2 F (l + Fo) || / || —L—^-w + + - using Q 

1 + (C + fc - As 2 )e 2 ^ 

< w. 

2s 

Since s is large, C + k — As 2 is negative so that using (|2.1|) we have 

I Cs+fcs-As 3 

C — A? iS — I — , . 

to < 3 - — w. (3.5) 



Now 



3 Cs + ks- As 3 _ 3C + Cs - 2ks 
+ C- ks+ fs 3 ~ C- A:s+ fs 3 

/9C , 3C 6Jfe\ 

s — 8 
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so that setting C"(A) := - A 6fc we obtain the estimate 



Thus (|3.5|) implies that 



w C"(A) , 

u; < 1 —s 6 w 

s 2 



which multiplied by s £ gives 

d , i e x i e Z' i C"( A ) s 
-r(s 1_e w < s 1-£ u; -es" 1 + — — s -3 
as \ 2 

< since s is large. 
Thus the function s t— * s 1 ~ £ w(s) is decreasing on This implies that 



to 



by assumption on t\ and so 



t ^ - 1+£ 



w(t) < w I -J ■ (3-7) 

This estimate holds only for characteristics which satisfy (|.3.2|) . but this is 
sufficient to conclude that 

P+(t)<w (J^\ ,te[h,t 2 ], 

in contradiction to the choice of t\. The estimate on P + is now established. 
The analogous arguments for characteristics with w < yield the assertion 
for P_. 

Next we consider characteristics which are not in the support of /. We can 
rewrite the inequality 1)3. 4|) for s £ [to,t] and u>(s) > : 

.2^, „ n n °, n e2M n2 ^ 1 1 + (A; - As 2 )e 2 ^ 

Prom 1)2.1)1 and 1)3.6)1 it follows that 1+ ^~ 2 ^ s ^ e M < 0. Using the estimate 
(|3~T|) on P + and ((277)1 we obtain 

w < C S 2£ " 3 — . 

2w 

Hence 

as 
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Integrating this over [to,t] yields 

w 2 (t) <C, t> t . (3.8) 

The analogous arguments for characteristics outside the support of / with 
w < yield the same estimate. Thus by (|3.7|) . (|3.1[) and (|3.<S[) we can state 

Proposition 3.1 For any characteristic (r, w, F), for any solution of Einstein- 
Vlasov system with positive cosmological constant and plane or hyperbolic 
symmetry written in areal coordinates and with initial data as in Proposi- 
tion^^ 

\w(t)\ <C, t> t , 
where the positive constant C depends on the initial data. 

3.2 Geodesic completeness 

Let ]t_,t + [9 t i — > (x a (t) , p@ (t)) be a geodesic whose existence interval is 
maximally extended and such that x°(tq) = i(ro) = to fo r some To e]r_, t+[. 
We want to show that for future-directed timelike and null geodesies r + = 
+oo. Consider first the case of a timelike geodesic, i.e., 

OafBP ? 13 = ~m 2 ; p° > 

with m > 0. Since ^ = p° > 0, the geodesic can be parametrized by 
the coordinate time t. With respect to coordinate time the geodesic exists 
on the interval [to,oo[ since on bounded t- intervals the Christoffel symbols 
are bounded and the right hand sides of the geodesic equations written 
in coordinate time are linearly bounded in p , p 2 , p 3 . Recall that along 
geodesies the variables t, r, p°, w := e p 1 , F := t A [(p 2 ) 2 + sin 2 #(j> 3 ) 2 ] 
satisfy the following system of differential equations : 

*1 = e - V ^ = -Xp°w - ^"VCp ) 2 , ^ = (3.9) 
dr dr dr 



ti =p o^ = _ A(p o )2 _ 2e- x n'p° W - e~^\w 2 - e~ 2 H^F. (3.10) 
dr dr 

Along the geodesic we define w and F as above. Then the relation between 
coordinate time and proper time along the geodesic is given by 



dt 



o 



— = p- = e Vy/m 2 + w 2 + F/t 2 , 
dr 

and to control this we need to control w as a function of coordinate time. 
By (J2.1|) we have the estimate 

e M > Ct~ l , t > to. 
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Combining this with the estimate on w in Proposition ^. ll vields the following 
along the geodesic : 

dr > Ct~ l 

dt ~ ^m 2 + w 2 + F/t 2 ~ Vm 2 + C + F' 

Since the integral of the right hand side over [io>oo[ diverges, r+ = +00 
as desired. In the case of a future-directed null geodesic, i.e. m = and 
P°( T o) > 0, p° is everywhere positive and the quantity F is again conserved. 
The argument can now be carried out exactly as in the timelike case, im- 
plying that r + = +00. We have proven : 

Theorem 3.2 Consider initial data with plane or hyperbolic symmetry for 
the Einstein- Vlasov system with positive cosmological constant. Suppose that 
the regularity properties required in the statement of Proposition^^ are sat- 
isfied. If the gradient of R is initially past-pointing then there is a corre- 
sponding Cauchy development which is future geodesically complete. 

3.3 Determination of the leading asymptotic behaviour 

In this subsection we determine the explicit leading behaviour of A, fj,, A, 
fi, //, and later on we compute the generalized Kasner exponents and prove 
that each of them tends to 1/3 as t tends to +00. 

Let us recall briefly some of the relevant notation. Let I be a set of real 
numbers and t* a real number or infinity. The asymptotic behaviour of a 
function g defined on / as t — * i* is to be described. It will be compared with 
a positive function h(t), typically a power of t. The notation g(t) = 0(h(t)) 
as t — > t* means that there is a neighbourhood U of t* such that there is 
a constant C with \g(t)\ < Ch(t) for all t belonging to both I and U. The 
notation g(t) = o(h(t)) as t — ► i* means that g{t)/h{t) tends to as t — > t*. 

Now Q1.4j) can be written in the form 

^(*e" 2 ") = At 2 - k-8irt 2 p. (3.11) 
Integrating this over [to,t] yields 

t e - 2 » = (t e- 2M * o) + kt - + ^t 3 -kt- f 87rs 2 pds. (3.12) 

3 3 J to 

By (|3.7jl we have the following, where C is a positive constant and e e]0, 1[ 

w < Cr l+e for t > to- 
Using the expression (|1.8|) for p, this implies that 

p < cr 5+3e 
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so that 

8vrf 2 p < Ct- 3+3£ . (3.13) 
Assuming e < 2/3 we obtain, using (|3.12|) 

i.e., 

e~^ = ^t 2 (l + 0(t- 2 )). 

It follows that 

^ = yft" 1 (1 + 0(t- 2 )) • (3.14) 

Now by (fPjl. and using (j3~Tij) and the fact that 8trtp = 0(t~ 2+£ ), we have 

• 1 /A . 2)) l + /ce 2 ^ 
A = -(At + 8ntp)e 2 ^ 



2 V r ' 2t 



r 2 (1 + 0{t -2^ (At + 0(t -5H*)) - - - - ^- 2 (1 + 0(t- 2 )) 



2A v ' " ' y 2t 2t 

and hence 



3 



\ = t- 1 (l + 0(t- 2 )) . (3.15) 
Integrating this over [to , i] yields 

A = lni [1 + O ((hit)- 1 )] . (3.16) 

Next, using ifTTIJl . $£T3§ and (|3~TH) we have 

fx = -r 1 (I + 0(t~ 2 )) (3.17) 

and integrating this over [to,t] yields 

// = -lni [l + 0((lni)~ X )] . (3.18) 

Now (j3.16|) implies that 

e x = 0(t), 
the expression p. 9(1 of j implies that 

bl < cr 4+2e 

and thus using equation ((1.5(1 we obtain 

^ = 0(t- 3 + 2£ ). (3.19) 

We can now compute the limiting values of the generalized Kasner exponents 
namely 

Kl(t,r) __ iA(t,r) K 2 2 (i,r) _ Kl(t,r) _ 1 



#(t,r) tA(t,r) + 2 K(t,r) K(t,r) t\(t,r) + 2' 

14 



where K(t, r) = K%(t, r) is the trace of the second fundamental form K^j of 
the metric. We refer to [Hj for the computation. Using (|3.15|) . we see that 
as t tends to +00, each of those quantities tends to 1/3 uniformly in r G 1. 
We have proved the following : 

Theorem 3.3 Let (/, A, fi) be a solution of the Einstein- Vlasov system with 
plane or hyperbolic symmetry and A > given in the expanding direction. 
Then the following properties hold at late times : \3.15\) . \3.1b]) . \'J. 
Unto . UTWtl . with e e]0,2/3[ and 

hm hm lim MM = I. 

K(t,r) t^+oo K(t,r) K(t,r) 3 

This theorem shows how the de Sitter solution acts as a model for the 
dynamics of the class of solutions considered in this paper. For if we set 
A = lni, // = — mi and k = the spacetime obtained is the de Sitter 
solution. Thus the leading terms in the asymptotic expansions of the metric 
components are exactly the quantities defined by the de Sitter spacetime. 
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